A0=10;
A1=5;
f=6000;
pi=3.14;
w0=f*pi;
t=-1:0.00001:1;
% 调制信号mes为A1*cos(0.001*w0*t)
mes=A1*cos(0.001*w0*t);
subplot(321);plot(t,mes)
title('调制信号');
Y1=fft(mes);
subplot(322);plot(abs(Y1));title('调制信号频谱')
axis([198000,202000,0,1000000]);


% 载波信号Uc为A0*cos(w0*t)
Uc=A0*cos(w0*t);
subplot(323);plot(t,Uc)
axis([0,0.01,-15,15]);
title('载波信号');
Y2=fft(Uc);             %对载波信号进行傅里叶变换 
subplot(324);plot(abs(Y2));title('载波信号频谱');  % abs绝对值 
axis([5800,6200,0,1000000]); 

% AM调制信号波形
t=-1:0.00001:1;  
A0=10;  %载波信号振幅
A1=5;   %调制信号振幅
A2=3;   %已调信号振幅
f=3000;  %载波信号频率
w0=2*f*pi;
m=0.15;    %调制度
mes=A1*cos(0.001*w0*t);   %消调制信号
Uam=A2*(1+m*mes).*cos((w0).*t); %AM 已调信号  
subplot(325);
plot(t,Uam);
grid on;  % 打开网格
title('AM调制信号波形');
subplot(326);
Y3=fft(Uam);              %  对AM已调信号进行傅里叶变换  
plot(abs(Y3)),grid;
title('AM调制信号频谱');
axis([5950,6050,0,500000]); 
